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ABSTRACT 

Glycolaldehyde is a key molecule in the formation of biologically relevant molecules such as ribose. We report its detection with 
the Plateau de Bure interferometer towards the Class 0 young stellar object NGC 1333 IRAS2A, which is only the second solar- 
type protostar for which this prebiotic molecule is detected. Local thermodynamic equilibrium analyses of glycolaldehyde, ethylene 
glycol (the reduced alcohol of glycolaldehyde) and methyl formate (the most abundant isomer of glycolaldehyde) were carried out. 
The relative abundance of ethylene glycol to glycolaldehyde is found to be ~5 - higher than in the Class 0 source IRAS 16293- 
2422 (— 1), but comparable to the lower limits derived in comets (> 3-6). The different ethylene glycol-to-glycolaldehyde ratios in 
the two protostars could be related to different CH 3 OH:CO compositions of the icy grain mantles. In particular, a more efficient 
hydrogenation on the grains in NGC 1333 IRAS2A would favor the formation of both methanol and ethylene glycol. In conclusion, it 
is possible that, like NGC 1333 IRAS2A, other low-mass protostars show high ethylene glycol-to-glycolaldehyde abundance ratios. 
The cometary ratios could consequently be inherited from earlier stages of star formation, if the young Sun experienced conditions 
similar to NGC 1333 IRAS2A. 

Key words, astrochemistry - astrobiology - stars: formation - stars: protostars - ISM: molecules - ISM: individual object 
(NGC 1333 IRAS2A) 


1. Introduction 


The inner regions of low-mass protostars are known to harbor a 
rich complex organic chemistry characterized by the presence of 
molecules such as methyl formate (CH3OCHO), dimethyl ether 
(CH3OCH3), and ethyl cyanide (C2H5CN) (e.g., Cazaux et al.| 


2003 1 |Bottinelli et al.|2004a[ |Bisschop et al.|2008 ” To differen¬ 

tiate them from the hot cores present in high-mass star-forming 
regions, they were called hot corinos ( |Ceccar elli 2004; Bott inelli] 
et al.|2004b| ). These complex organic molecules are thought to be 
efficiently formed on grains and then released into the gas phase 
in the hot corino by thermal desorption (e.g., jGarrod et al.j2008| 
Herbst & van Dishoeckj 2009[ ). Some of these complex organic 
molecules are particularly interesting because of their supposed 
role in the emergence of life. Indeed, the detection of so called 
prebiotic molecules in low-mass star-forming regions indicates 
that they can form early during the star formation process and 
thereby be available for possible later incorporation into solar 
system bodies, e.g., comets. 

Glycolaldehyde (CHiOHCHO) is one of these prebiotic 
molecules: it is a simple sugar-like molecule and under Earth¬ 
like conditions the first product in the formose reaction leading 
to the formation of ribose, an essential constituent of ribonucleic 
acid (RNA) (e.g., |Zubay & Mui|2001|rJalbout et al.|2007] l. Gly¬ 


colaldehyde was first detected towards the Galactic center (Sgr 
B2(N): |Hollis et al.|2000l [2001][2004] |Halfen et al.f2006l [Bet] 


* Based on observations carried out with the IRAM Plateau de Bure 

[Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 


loche et al. 2013 molecular clouds: Requena-Torres et al. 2008). 


Later it was shown to be present in the high mass star-forming re¬ 
gion G31.41+0.31 ( Beltran et al.|2009), in the in termediate mass 
protostar NGC 7129 FIRS 2 ( Fuente et al.|2014] l, and even in the 
hot corinos of the Class 0 protostellar binary, IRAS 16293-2422 
(hereafter IRAS 16293, [Jprgensen et al.|2012j l. This indicates that 
this molecule can be synthesized relatively early in the environ¬ 
ments of solar-type protostars. Furthermore, glycolaldehyde can 
easily survive during impact delivery to planetary bodies, and 
impacts can even facilitate the formation of even more complex 
molecules ( |McCaffrey et al.|2014) >. 


Similarly to other complex organic molecules, the formation 
of glycolaldehyde is thought to occur on grains. In particular, a 
gas phase formation was excluded by Woods et al. (2 012||20l3] l, 
as the produced abundances are too low compared with the ob¬ 
servations. Several grain surface formation pathways were pro¬ 
posed in the literature. | Woods et al.| ( j20T2] i modeled their effi¬ 
ciency and showed that the formation by the reaction CH3OH 
+ HCO would be very efficient, but that, from chemical consid¬ 
erations, H3CO + HCO could be more feasible. Another proba¬ 
bly efficient way to form glycolaldehyde would be through HCO 
dimerization (HCO + HCO —> HOCCOH) followed by two suc¬ 
cessive hydrogenations ( (Woods et al.||2013) . A recent experi¬ 
mental study based on surface hydrogenations of CO seems to 
confirm this pathway (Fedoseev et al. submitted I. 


A related species to this prebiotic molecule is ethylene gly¬ 
col ((CH 2 OH) 2 ). More commonly known as antifreeze, it is the 
reduced alcohol of glycolaldehyde. This molecule was tenta¬ 
tively detected towards IRAS 16293 with one line of the gGg' 
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conformer (J 0 rgensen et al. 


2012). Interestingly, the aGg' con- 


former of ethylene glycol (the conformer of lowest energy) was 
detected in three comets, Hale-Bopp, Lemmon, and Lovejoy, 
while glycolaldehyde was not, leading to a lower limit of 3-6 
for the (CH 20 H) 2 /CH 20 HCH 0 abundance ratio ( |Crovisier et al.| 
2004 1 Biver et al.|2014 1. Ethylene glycol was also detected in the 


Murchison and Murray carbonaceous meteorites, while the pres¬ 


ence of aldehyde sugars have not been reported yet (Cooper et 

|al .|2()0 1 


NGC 1333 IRAS2A (hereafter IRAS2A) is another of these 
famous hot corinos. In particular, methyl formate, the most abun¬ 
dant isomer of glycolaldehyde, was detected towards this source 
by |j 0 rgensen et al.| ( |2005a| ) and |Bottinelli et al.| ( |2007| i. More 
recently, ethylene glycol was detected in the framework of the 
CALYPSO program carried out with the IRAM Plateau de Bure 
Interferometer (PdBI) by [Maury et al.j ( [2014] ). We here report 
the detection of glycolaldehyde towards the same low-mass pro¬ 
tostar, and present an analysis of the relative abundances of these 
three species. 


2. Observations 


This work is based on several separate programs carrying out 
observations of the solar-type protostar IRAS2A with the PdBI. 
Four spectral ranges (84.9-88.5, 223.5-227.1, 240.2-243.8, 
and 315.5-319.1 GHz) were covered with the WIDEX corre¬ 
lator at a spectral resolution of 1.95 MHz (dv = 6.8 km s 1 at 86 
GHz, dv = 2.6 km s _1 at 225 GHz, dv = 2.4 km s _1 at 242 GHz, dv 
= 1.8 k m s 1 at 317 GHz) and reduced with the GILDAS0 soft¬ 
ware. The synthetized beam sizes obtained with natural weight¬ 
ing are about 3.0" x3.0" at 86 GHz, 1.2" x 1.0" at 225 GHz, 
1.4" x 1.0" at 242 GHz, and 0.9" x 0.8" at 317 GHz. The dust 
continuum fluxes at 0.9 and 1.3 mm are consistent with previ¬ 
ous measurements (e.g., J 0 r gensen et al.|[2007[ |Persson et al.| 
2012). The absolute calibration uncertainty for each dataset is 


about 20%. Additional information about the observations and 
their reduction can be found in jCoutens et al.| ( |2014| ) and |Persson| 
et al. (2014[). The 3mm data are from Wampfler (priv. comm.). 


Using the CASSIS^] software, we detected 8 lines of glyco¬ 
laldehyde, 31 lines of the aGg' conformer of ethylene glycol, 
and 26 lines of methyl formate (see Table[Tji. The glycolaldehyde 
and methyl formate transitions are taken from the JPL spectro¬ 
scopic database ( |Pickett et al.||1998) , while the ethylene glycol 
transitions are from the CDMS catalogue (Muller et al. 2001 [ 
2005 1). The predictions are based on experimental data from|But 


ler et al.| ( |2001| ), |Widicus Weaver et aT7|(]2005[> a nd |Carroll et al. 
(|2010|) for glycolaldehyde, Christen et al. ( 1995|) and Christen & 


Muller] < [2003| ) for ethylene glycol, and |Ilyushin et al.| ( |2009| ) for 
methyl formate. The frequencies of five of the detected glyco¬ 


laldehyde lines were directly measured in the laboratory (Butler 
et al.] |2001| ). Some of the lines result from a blending of sev¬ 
eral transitions of the same species. The lines that are strongly 
blended with other species are not listed in Table[l] All three 
species are emitted very compactly at the position of the contin¬ 
uum peak (o' 2 ooo- 03 h 28 m 55 ? 57 , (>2000 = 31°14'37'.'l). The angular 
sizes obtained with a circular Gaussian fit in the (u,v) plane vary 
from a point source to a maximum of 1 " depending on the tran¬ 
sition. The line fluxes listed in Table|]]were measured at the con¬ 
tinuum peak position with the CASSIS software using a Gaus¬ 
sian fitting method (Levenberg-Marquardt algorithm). The lines 
that are contaminated in the wings by other transitions are con- 


1 http://www.iram.fr/IRAMFR/GILDAS/ 

2 http://cassis.irap.omp.eu 
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Fig. 1. Rotational diagrams for glycolaldehyde, ethylene glycol, and 
methyl formate. 


sequently fitted with a sum of Gaussians. We carefully checked 
that the derived full widths at half maximum ( FWHM ) are con¬ 
sistent with the other line measurements. The average FWHM 
is about 4.5 kms -1 at 317 GHz, and 5.0 kms -1 at 225 and 242 
GHz. The widths of the methyl formate lines at 87 GHz are quite 
broad (-12 km s“ 1 ). It is consequently difficult to completely ex¬ 
clude an extra flux contribution from other species. The variation 
of FWHM with the frequency can be explained by the spectral 
resolution of the observations that decreases towards the lower 
frequencies. 


3. Results 

We carried out a local thermodynamic equilibrium analysis of 


the three species through the rotational diagram method (Gold¬ 
smith & Langer 1999|. We consider that the lines are emitted 


in a region of 0.5" size, which is the average size derived for 
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Fig. 2. Observed lines of CH 2 OHCHO towards the protostar NGC 1333 IRAS2A (in black). The best-fit model for CH 2 OHCHO (see Table [ 2 } 
is shown in red solid lines. The line in the third upper panel is blended with an unidentified species. The contribution of the CH 3 OCHO lines is 
indicated with green lines, and the model including both CTFOHCHO and CH 3 OCHO can be seen in red dotted lines. 


the methyl formate lines when fitting circular gaussians in the 
( u,v ) plane (see also |Maury et al.|2014[). It is also similar to what 
we found for deuterated water (Coutens et a l.|2014) . It also cor¬ 
responds to the expected size of the region where the tempera¬ 
ture increases above ~ 100 K according to dust radiation trans¬ 


fer models of the envelope (Jprgensen et al. 2002 1 and where 
the complex molecules and deuterated water should sublimate 
from the grains. The line fluxes that result from a combination 
of several transitions of the same species are used in the rota¬ 
tional diagrams, unless the transitions have different E up val¬ 
ues. For glycolaldehyde, we include two lines slightly blended 
with some methyl formate transitions after subtraction of the 
predicted flux contribution from methyl formate. As the best- 
fit model for methyl formate reproduces extremely well the ob¬ 
servations, the final fluxes of the glycolaldehyde lines can be 
trusted, which is also confirmed by their alignment with the other 
points in the rotational diagram of glycolaldehyde (see upper 
panel in Figure [I). Assuming a source size of 0.5", we derive 
column densities (with lcr uncertainties) of 2.4^ ® x 10 15 cm~ 2 , 

and excitation tem- 
7 K for glycolalde- 


i+0.1 
-0.1' 


316, 


2+0.3 

-0.3 


x 10 16 cm 


> + ®K, 133+?! K, and 126 +9 


peratures of 103! 20 ii., 14 j 

hyde, ethylene glycol, and methyl formate, respectively. Within 
the uncertainty range, the excitation temperature seems to be 
similar between the three species (~ 130 K), which is consistent 
if the three species arise from a same region. We checked, for 
each species, that there is no line flux overpredicted by the model 
anywhere in the four datasets. In the case of ethylene glycol, the 
model shows an overproduced flux for some transitions, espe¬ 
cially the lines (240.778, 241.545, 241.860, and 316.444 GHz) 
that correspond to the four lower points in the rotational diagram 
(see middle panel in Figure[T]i. A model with a column density of 
1.1 x 10 16 cm ~ 2 would be sufficient to produce line fluxes consis¬ 
tent with these observations. Table [2] summarizes the parameters 
used for the line modeling of the three species that can be seen 
in Figures [2] [3j and [4] According to these models, all lines are 
optically thin (r < 0 . 1 ). 


Although no other species than glycolaldehyde is found at a 
frequency of 240890.5 MHz, the line (E up = 52 K) is probably 


Table 2. Parameters used to compute the synthetic spectra of glyco¬ 
laldehyde, ethylene glycol, and methyl formate. 


Molecule 

Source 
size (") 

7?x 

(K) 

N 

(cm -2 ) 

t+SR 

(kms -i ) 

ch 2 ohcho 

0.5 

130 

2.4 x 10 1S 

7.0 

aGg'-(CH 2 OH ) 2 

0.5 

130 

1.1 xlO 16 

7.0 

CH 3 OCHO 

0.5 

130 

4.8 xlO 16 

7.0 


Notes. The FWHM used for the line modeling are 4.5, 5.0, 5.0, and 7.0 
km s _1 for the data at 317, 242, 225, and 86 GHz, respectively. 


blended with an unidentified species: the predicted flux is com¬ 
pletely underproduced with respect to the observations, and it 
cannot be due to a different excitation in the cold gas, as a line of 
glycolaldehyde at 243232.21 MHz (£ up = 47 K) - blended with 
a bright CH 2 DOH line in the red-shifted part of the spectrum 
and also potentially blended with a DCOOH line (E up = 106 K, 
Aq = 1.35 x 10 ~ 4 s -1 ) - would have a higher flux inconsistent 
with the observed one. 


4. Discussion 

The relative abundances of the three species are derived 
from the column densities in Table [2] and compared with 
other star-forming regions and comets in Table [3] The 
(CH 2 OH) 2 /CH 2 OHCHO abundance ratio of -0.3-0.5 previ- 


ously derived in IRAS 16293 by Jdrgensen et al. ( 

2012 ) was 

revised. Indeed, the assignment in Jprgensen et al. 

2012 ) was 

based on only one line of the gGg' conformer of ethy 

ene glycol 

about 200 cm 1 (~290 K, 

Muller & Christen [2004 

above the 


lowest-energy aGg' conformer - and thus tentative. An analy¬ 
sis from observations of 6 transitions of the lower energy con- 
former from ALMA Cycle 1 observations at 3 mm (4 spec¬ 
tral windows at 89.48-89.73, 92.77-93.03, 102.48-102.73 and 
103.18-103.42 GHz; Jprgensen et al. in prep.) results in a higher 
ethylene glycol-to-glycolaldehyde abundance ratio of 1.0+ 0.3. 
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Table 3. (CH 2 OH) 2 /CH 2 OHCHO, CH 3 OCHO/CH 2 OHCHO, and CH 3 OCHO/(CH 2 OH) 2 column density ratios determined in different objects. 


Source 

f C H 2 O H) 2 /C H 2 O H C H O" 

CH 3 OCHO/CH 2 OHCHO C H 3 OC H 0/ f C H 2 O H ) 2 a 

References 



Class 0 protostars 


NGC 1333 IRAS2A 

~5 


~20 ~4 

1 

IRAS 16293-2422 

~1 


~13 ~13 

2 




Comets 


C/1995 Ol (Hale-Bopp) 

>6 


>2 ~ 0.3 

3 

C/2012 F6 (Lemmon) 

> 3 


<0.7 

4 

C/2013 R1 (Lovejoy) 

>5 


<0.6 

4 


High- and intermediate-mass star forming regions 


Sgr B2(N) 

0.7-2.2 b 


~52 c ~30 

5,6,7 

G34.41+0.31 



<34 

8 

NGC7129 FIRS2 

~2 


~40 ~20 

9 


Molecular clouds in the Central Molecular Zone 


G-0.02, G-0.11, and G+0.693 

— 1.2—1.6 


~3.3-5.2 ~2.5-4.3 

10 

Notes. (a) (CH 2 OH) 2 refers to the aGg' conformer. <fc) Large scale emission. Hot core emission. 

References. 1) this study, 2) Jprgensen et al. (2012, in prep.). 3)|Crovisier et al. 

(2004), 4) Biver et al. (2014), 5) Hollis et al. (2001), 6) Hollis eh 


This new estimate is consistent with the ratio expected be¬ 
tween the aGg' and gGg' conformers under thermal equilib¬ 
rium conditions at 300 K, the excitation temperature of glyco- 
laldehyde derived in IRAS 16293 ( |j0rgensen et al.||2012| >. The 
(CH20H)2/CH20HCH0 abundance ratio in IRAS2A is esti¬ 
mated at 5.5 + 1.0 if we consider the column densities derived 
from the rotational diagrams. It is however slightly lower (4.6) 
if we use the column density of ethylene glycol of 1.1 x 10 16 
cm~ 2 that does not overproduce the peak intensities of a few 
lines (see Fig. [3j. The (CfTOH^/CIFOHCHO abundance ratio 
is consequently a factor ~5 higher than in the Class 0 protostar 
IRAS 16293. It is also higher than in the other star-forming re¬ 
gions (see Table [3J, but comparable to the lower limits derived 
in comets (> 3-6). This indicates that the glycolaldehyde chem¬ 
istry may vary among hot corinos in general. It is possible that, 
like IRAS2A, other very young low-mass protostars show high 
(CH20H)2/CH20HCH0 abundance ratios, in agreement with 
the cometary values. The CH3OCHO/CH2OHCHO column den¬ 
sity ratio found in IRAS2A (~ 20) ranges between the values 
derived in the molecular clouds from the Galactic Center (~ 3.3- 
5.2) and the high-mass star-forming regions (~ 40-52). A lower 
limit of 2 was derived for the Hale-Bopp comet. 

In contrast to IRAS 16293, the (CH 2 OH) 2 /CH 2 OHCHO 
abundance ratio in IRAS2A is comparable to the lower lim¬ 
its in comets. To explain these different abundance ratios in 
IRAS2A and IRAS 16293, two scenarios are possible: either the 
(CH20H)2/CH20HCH0 ratio is similar in the grain mantles of 
low-mass protostars and it evolves in the gas phase after the sub¬ 
limation of the molecules in the hot corinos, or this ratio was 
already different in the grain mantles of the two protostars. 


In the first scenario, if we assume that the 
(CthOH^/ChkOHCHO increases until it reaches the cometary 
value, it would mean that glycolaldehyde can easily be de¬ 
stroyed in the gas phase of the warm inner regions. Another 
possibility would be that ethylene glycol can form efficiently 
in the gas phase, but complex organic molecules are generally 
difficult to form with high abundances in the gas phase. If the 
evaporation temperature of ethylene glycol is higher than glyco¬ 
laldehyde, as assumed in the chemical model of Garrod ( 2013| >, 
ethylene glycol would desorb later than glycolaldehyde and the 
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(CfFOH^/CIFOHCHO abundance ratio would consequently 
increase with time (until the two molecules have completely 
desorbed). This chemical model however predicts an abundance 
of glycolaldehyde significantly higher than those of ethylene 
glycol and methyl formate, which is inconsistent with the ratios 
derived in IRAS2A. More theoretical and experimental work 
would be needed to make the case that these hypotheses are 
plausible. 


In contrast, experimental studies based on irradiation of ices 
show that the second scenario is likely. Such studies show that 
glycolaldehyde, ethylene glycol, and methyl formate can be syn¬ 
thesized by irradiation of pure or mixed meth anol (CH3OH) ices 
(Hudson & Moore 2000, Oberg et al. 2009). Interestingly, the 
(CH20H)2/CH20HCH0 abundance ratio is found to be depen¬ 
dent on the initial ice composition as well as the ice temper¬ 
ature during the UV irradiation. The CH30H:CO ratio in the 
ices is a key parameter: for irradiated 20 K ices a composi¬ 
tion of pure CH3OH leads to a (CH20H) 2 /CH20HCH0 ratio 
higher than 10, while a CH30H:CO 1:10 ice mixture produces 


a (CH20H) 2 /CH20HCH0 ratio lower than 0.25 (Oberg et al. 
12009 ). The difference found between IRAS 16293 and IRAS2A 
could then be related to a different grain mantle composition 
in the two sources. If the CH30H:CO ratio in the grain man¬ 
tles of IRAS 2A was higher than in IRAS 16293, a higher 
(CH20H) 2 /CH20HCH0 abundance ratio would be expected 
according to the laboratory results. In fact, the CH3OH gas- 


phase abundance in the inner envelo 

pe is found to be higher in 

IRAS2A(~4x lO~ 7 ,|j0r2ensen etal. 

2005b (than in IRAS 16293 

(~ 1 x 10- 7 , Schoier et al. 2002), whi 

e the CO abundance is rel- 

atively similar (~(2-3)xl0 , J0rgensen et al. 2002 Schoier 


|et al.||2002) . This could consequently be the result of the des 
orption of ices with a higher CH30H:CO ratio in IRAS2A than 
IRAS 16293. The question then arises: how can CH3OH be more 
efficiently produced on grains in IRAS2A than in IRAS 16293? 
Several scenarios are possible: i) The initial conditions may play 
an important role in the CH 3 OH:CO ratio. In particular, exper¬ 
iments and simulations show that the efficiency of CH3OH for¬ 
mation through CO hydrogenation on the grains is dependent 
on temperature, ice composition (CO:H20), and time (|Watan- 


|abe et alTj|2004[ [Fuchs et ai.||2009] l. ii) The coiiapse timescaie 
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radius (All) 

Fig. 5. Black: Ho density (solid line) and temperature (dotted line) pro¬ 
files of the protostar IRAS2A from |j0rgensen et aTT| i }2002[ >. Red: H 2 den¬ 
sity (solid line: power-law model, dashed line: Shu-like model) and tem¬ 
perature (dotted line) profiles of the protostar IRAS16293 from |Crimier| 
[etaL| ( [20T0| . 


was longer in IRAS2A than in IRAS 16293, enabling to form 
more CH3OH. iii) The H 2 density in the prestellar envelope of 
IRAS2A was lower than that of IRAS 16293. Indeed a less dense 
environment would lead to a higher atomic H density and con¬ 
sequently to a higher efficiency of CO hydrogenation. This was 
proposed by Maret et al. ( 2004| > and Bottinelli et al.| ( 2007| l to 
explain an anti-correlation found between the inner abundances 
of H 2 CO and CH3OH and the submillimeter luminosity to bolo- 
metric luminosity (L smrn /L ho |) ratios of different low-mass pro¬ 
tostars. The L smm /Lboi parameter is interpreted as an indication 
of different initial conditions, rather than an evolutionary param¬ 
eter in this context (Maret et al.]|2004j>. The L smm /Lboi ratios of 
IRAS2A (~ 0.00 5, |Karska et al.|2013| > and IRAS 16293 (-0.019, 
Froebrich|2005 1 are consistent with this hypothesis. The current 


H2 density profiles of these two sources are also in agreement 
with this scenario, if they keep the memory of the prestellar con¬ 
ditions. The density derived in the outer envelope of IRAS2A 
with a power-law model (Jprgensen et al. 20021 is lower than the 
density derived in IRAS 16293 by Crimier et al. ( |2010j > whether 
it be for a Shu-like model or a power-law model, while the tem¬ 
perature profiles are relatively similar (see Figure [5]). Along the 
same lines,|Hudson et al. (2005[ i showed with proton irradiation 
experiments that glycolaldehyde is more sensitive to radiation 
damage than ethylene glycol. Irradiation would be more impor¬ 
tant in less dense envelopes, which would also be consistent with 
a less dense prestellar envelope in IRAS2A. A recent experiment 
by Fedoseev et al. (submitted) shows that these two species can 
also be synthesized by surface hydrogenations of CO molecules 
in dense molecular cloud conditions. They do not directly form 
from CH3OH, but the results of this experiment show that sim¬ 
ilarly to CH3OH that results from successive hydrogenations of 
CO, ethylene glycol forms by two successive hydrogenations of 
glycolaldehyde. This is consequently in agreement with the pro¬ 
posed scenario. 

In conclusion, the (CH 2 OH) 2 /CH 2 OHCHO abundance ra¬ 
tio measured in low-mass protostars can be different from one 
source to another, and possibly consistent with cometary values. 
In some cases, the (CH 2 OH) 2 /CH 2 OHCHO ratios determined in 
comets could consequently be inherited from early stages of star 


formation. Such a difference between low-mass protostars could 
be related to a different CH30H:CO ratio in the grain mantles. 
A more efficient hydrogenation (due for example to a lower den¬ 
sity) on the grains would lead to higher abundances of CH3OH 
and (CH 2 OH) 2 . A determination of (CH 2 OH) 2 /CH 2 OHCHO ra¬ 
tios in larger samples of star-forming regions could help under¬ 
stand how the initial conditions (density, molecular cloud, ...) 
affect their relative abundances. 
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Fig. 3. Observed lines of aGg'-(CH 2 OH ) 2 towards the protostar NGC 1333 IRAS2A (in black). The best-fit model for aGg'-(CH 2 OH ) 2 (see Table 
[2) is shown in red solid lines. The contributions of the CH 3 OCHO and CH 2 OHCHO lines are indicated with green and blue lines, respectively. 
The model including aGg'-(CH 2 OH) 2 , CH 2 OHCHO, and CH 3 OCHO can be seen in red dotted lines. 
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Fig. 4. Observed lines of CH3OCHO towards the protostar NGC 1333 IRAS2A (in black). The best-fit model for CH3OCHO (see Table[2| 

is shown in red lines. 
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Table 1 . CH 2 OHCHO, aGg'-(CH 2 OH) 2 and CH 3 OCHO transitions observed towards NGC 1333 IRAS2A. 


Species 

Transition 

Frequency 

(MHz) 

£u P 

(K) 

Aij 

(s- 1 ) 

£up 

Flux 

(Jy km s -1 ) 

RD“ 

CH 2 OHCHO 

13io,3 — 13g,4 (T=0) 

240366.34* 

111.3 

1.2 x 10~ 4 

27 

0.055 

Y 


1310 , 4 - 139,5 ( v-0 ) 

240366.34* 

111.3 

1.2 xlO- 4 

27 



CH 2 OHCHO 

12io, 2 -129,3 (u= 0 ) 

240482.78* 

104.0 

1.0 xlO" 4 

25 

0.038 

Y 


12 io ,3 —129,4 (v-0) 

240482.78* 

104.0 

1.0 xlO- 4 

25 



CH 2 OHCHO 

115 . 6 - 104,7 (u= 0 ) 

240890.46 

51.9 

1.8 x 10~ 4 

23 

0.663 

N b 

CH 2 OHCHO 

22 2 ,20-2l3,i9 ( v-0 ) 

241131.84 

142.8 

2.8 x 10~ 4 

45 

0.061 

Y 

CH 2 OHCHO 

23 2 , 22 — 22 1>2 i ( v - 0 ) 

242239.09 

146.2 

3.5 x 10~ 4 

47 

0.128 

Y 

CH 2 OHCHO 

240 , 24 - 23!,23 (v=0) 

242957.72* 

148.2 

4.2 x 10~ 4 

49 

0.258 

Y 


24i.24-23o.23 (v=0) 

242957.98* 

148.2 

4.2 x 10" 4 

49 



CH 2 OHCHO 

19i3,7 — 19i 2 ,8 (tf=0) 

315941.48* 

208.1 

3.2 x 10- 4 

39 

0.074 

Y 


19i3,6 — 19i 2 ,7 ( v - 0 ) 

315941.48* 

208.1 

3.2 x 10- 4 

39 



CH 2 OHCHO 

ll8,4-10 7 ,3 («=0) 

317013.88* 

75.5 

6.6 x 10- 4 

23 

0.363 

Y r 


118,3 — 107,4 (v=0) 

317013.90* 

75.5 

6.6 x 10- 4 

23 



CH 2 OHCHO 

275 , 23 - 264,22 (v-0) 

317850.44 

226.2 

4.5 x 10~ 4 

55 

0.117 

Y c 

aGg'-(CH 2 OH ) 2 

216,16 (f=l)—206,15 (f=0) 

223741.66 

132.0 

2.5 x 10“ 4 

387 

0.212 

Y 

aGg'-(CH 2 OH ) 2 

216,15 (y=l)-206,14 («=0) 

224405.85 

132.1 

2.5 x 10~ 4 

301 

0.197 

Y 

aGg'-(CH 2 OH ) 2 

24o,24 (y=l)-23i,23 (v=0) 

224511.70* 

136.8 

5.4 x 10~ 5 

441 

0.145 

N 


24i,24 G>=1)-23o,23 (^=0) 

224512.74* 

136.8 

5.4 xlO' 5 

343 



aGg'-(CH 2 OH ) 2 

213,18 (t ,= l)— 2 O 337 (u= 0 ) 

225688.94 

121.3 

2.4 xlO" 4 

301 

0.195 

Y 

aGg'-(CH 2 OH ) 2 

223,20 (tt— 1 )— 2 I 319 (u= 0 ) 

225929.69 

127.8 

2.5 x 10- 4 

315 

0.349 

N f 

aGg'-(CH 2 OH ) 2 

22 5 ,i7 (o= 0 )- 21 5 ,i 6 (it=l) 

226095.96 

138.2 

2.6 x 10~ 4 

315 

0.303 

N 

aGg'-(CH 2 OH ) 2 

22 2 , 2 o (o=l)- 21 2 ,i 9 (u= 0 ) 

226561.99 

127.7 

3.0 xlO' 4 

405 

0.329 

Y 

aGg'-(CH 2 OH ) 2 

25i,25 (u=0)-24i, 24 (tt=l) 

226643.30* 

147.7 

2.9 x 10~ 4 

357 

0.470 

Y 


25o,25 (h=0)—24o, 2 4 (u=1) 

226643.46* 

147.7 

2.9 x 10- 4 

459 



aGg'-(CH 2 OH ) 2 

25i, 25 (!t=l)-24i, 24 ( v = 0 ) 

240778.12* 

148.0 

3.4 x 10~ 4 

459 

0.360 

Y 


25o,25 (u=1)-24o, 24 (f=0) 

240778.30* 

148.0 

3.4 x 10~ 4 

357 



aGg'-(CH 2 OH ) 2 

24g,i7 (ti=0)-238,i6 (u=l) 

240807.88 

179.2 

3.0 xlO" 4 

441 

0.340 

N e 

aGg'-(CH 2 OH ) 2 

24g,i6 (ti=0)-23g,i5 (f — 1) 

240828.89 

179.2 

3.0 xlO" 4 

343 

0.149 

Y 

aGg'-(CH 2 OH ) 2 

245, 2 o (O=0)-23 5 ,19 (v=l) 

241291.27 

160.7 

3.1 xlO- 4 

441 

0.196 

Y 

aGg'-(CH 2 OH ) 2 

24 7 jg (y=0)-23 7 ,i7 (u=l) 

241545.26 

172.1 

3.1 xlO- 4 

441 

0.151 

Y 

aGg'-(CH 2 OH ) 2 

246,19 (y=0)-236,ig (i>=l) 

241860.73 

166.0 

2.8 x 10- 4 

441 

0.137 

Y 

aGg'-(CH 2 OH ) 2 

23is,g (p=1)-22i 5 ,7 (o=0) 

242244.69* 

246.4 

2.0 xlO" 4 

329 

0.474 

N 


23i5,9 (y=l)- 22 i 5 ,g (u= 0 ) 

242244.69* 

246.4 

2.0 xlO" 4 

423 




236,17 (!t=l)- 22 6 ,n (i>=l) 

242245.62* 

154.6 

1.1 x 10~ 5 

329 




2314,9 (y=l)- 22 i 4 ,g (u=0) 

242246.34* 

232.2 

2.2 x 10- 4 

329 




2314,10 (tt=l)-22i4, 9 (v-0) 

242246.34* 

232.2 

2.2 x 10- 4 

423 



aGg'-(CH 2 OH) 2 

23i3,io (u=1)-22i3,9 (d=0) 

242277.72* 

218.9 

2.4 xlO- 4 

329 

0.275 

N 


23i3.ii (^=1)-22i3,io (v-0) 

242277.72* 

218.9 

2.4 xlO' 4 

423 



aGg'-(CH 2 OH) 2 

23io,i4 (u=l)-22io,i3 (u=0) 

242656.22* 

185.2 

2.8 x 10- 4 

423 

0.305 

Y 


2310.13 (y=l)-22io,i 2 (v-0) 

242656.24* 

185.2 

2.8 x 10- 4 

329 



aGg'-(CH 2 OH) 2 

239,15 ( 11 = 1 )— 229,14 (u=0) 

242947.99* 

175.9 

3.0 xlO" 4 

423 

0.291 

Y 


239,14 ( v = 1 )- 22 913 (u=0) 

242948.59* 

175.9 

3.0 xlO" 4 

329 



aGg'-(CH 2 OH) 2 

235,19 (n=l)- 22 5 ,ig (u=0) 

243636.57* 

149.1 

3.4 xlO" 4 

423 

0.287 

Y 

aGg'-(CH 2 OH) 2 

318,23 (ll—0)—308,22 (n=l) 

315671.33 

276.6 

7.0 xlO" 4 

567 

0.314 

N f 

aGg'-(CH 2 OH) 2 

317,25 (n—0)— 307,24 (u=l) 

315892.11 

269.6 

6.8 x 10- 4 

441 

0.144 

N d 

aGg'-(CH 2 OH) 2 

309,21 (u=l >-299,20 (u-0) 

315961.89 

269.4 

7.0 xlO" 4 

549 

0.259 

Y 

aGg'-(CH 2 OH) 2 

34i,34 (u=0)-33o,33 (v=0) 

316444.07* 

268.5 

1.6 xlO' 4 

621 

0.078 

Y 


34o, 34 (u=0)—33 1,33 (u=0) 

316444.07* 

268.5 

1.6 xlO- 4 

483 



aGg'-(CH 2 OH) 2 

206,14 (tl=0)-19 5 ,15 (v=0) 

316698.08 

121.3 

5.1 x 10~ 5 

287 

0.029 

Y 

aGg'-(CH 2 OH) 2 

307,24 (tl=l)-29 7 ,23 (w=0) 

316868.23* 

254.4 

6.9 x 10- 4 

427 

0.445 

N r 


2016,4 (tl=0)-20i5,5 (w=0) 

316870.92* 

228.9 

3.8 xlO' 5 

287 




20i6,5 (y=0)-20i5,6 (v—0) 

316870.92* 

228.9 

3.8 xlO' 5 

369 



aGg'-(CH 2 OH) 2 

308,23 (n=l)-29 8 ,22 (u=0) 

316917.19 

261.4 

7.2 xlO' 4 

427 

0.299 

Y 

aGg'-(CH 2 OH) 2 

168,9 (ll=l)-15 7 ,8 (f=l) 

317054.30* 

98.6 

8.0 xlO" 5 

231 

0.116 

N 


168,8 (u=l)-15 7 ,9 (u=1) 

317055.36* 

98.6 

8.0 xlO" 5 

297 



aGg'-(CH 2 OH) 2 

2I4J8 (w=l)-204,i8(u=0) 

317267.77* 

122.1 

1.1 xlO- 5 

387 

0.312 

N 
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149.5 (b= 1)-13 8 ,6 (tt=l) 

317267.91* 

91.4 

1.1 xlO -4 

261 




149.6 (B=lH3y (w=l) 

317267.91* 

91.4 

1.1 xlO -4 

203 




149.5 (d=0)-13 8 ,6 (u=0) 

317268.88* 

91.4 

1.1 xlO -4 

203 




149,6 (u=0)-13 8 ,5 (p=0) 

317268.88* 

91.4 

1.1 xlO- 4 

261 



aGg'-(CH 2 OH) 2 

323,30 Crt— 1) 31 3,29 (v=0) 

317962.58 

257.2 

7.8 xlO- 4 

455 

0.310 

Y 

aGg'-(CH 2 OH) 2 

324.28 (lt=0)-314,27 (y=l) 

317982.56 

272.2 

4.8 x 10- 4 

455 

0.160 

Y 

aGg'-(CH 2 OH) 2 

35o,35 (y=0)-34o,34 0>=1) 

318433.40* 

284.1 

8.0 xlO- 4 

639 

0.639 

N 


351,35 (lt=0)-34i,34 (v= 1 ) 

318433.40* 

284.1 

8.0 xlO" 4 

497 



CH 3 OCHO 

73 4—63 3 E (u-f=0) 

87143.28* 

22.6 

7.7 x 10~ b 

30 

0.049 

N 


215,16-214,17 E (l>7’ = l) 

87143.65* 

342.0 

1.4 xlO " 6 

86 



CH 3 OCHO 

8o,8-2i,7 E (v t =1) 

87160.84* 

207.0 

1.3 xlO " 6 

34 

0.038 

N 


73 4—63 3 A (i>7=0) 

87161.28* 

22.6 

7.8 xlO " 6 

30 



CH 3 OCHO 

186,13-176,12 E (vr=0) 

224021.87* 

125.3 

1.5 xlO - 4 

74 

0.752 

Y 


186,13-176,12 A (Vt-0) 

224024.10* 

125.3 

1.5 xlO " 4 

74 



CH 3 OCHO 

I 85 . 14 -i 75.13 E (vt-0) 

224313.15 

118.3 

1.6 xlO " 4 

74 

0.351 

Y 

CH 3 OCHO 

193 , 17 - 183,16 E (vt-1) 

224491.31 

303.2 

1.7 xlO " 4 

78 

0.212 

Y 

CH 3 OCHO 

186,12-186,11 E (v T = 0) 

224582.35 

125.4 

1.5 xlO ' 4 

74 

0.463 

Y 

CH 3 OCHO 

186,12-176,11 A (vt-0) 

224609.38 

125.4 

1.5 xlO " 4 

74 

0.466 

Y 

CH 3 OCHO 

20 2 ,19—192.18 A (vt- 1) 

225372.22 

307.3 

1.7 xlO " 4 

82 

0.186 

Y 

CH 3 OCHO 

I 93 . 17 -I 83 .i 6 E (vt=0) 

225608.82 

116.7 

1.7 xlO " 4 

78 

0.426 

Y 

CH 3 OCHO 

I 85 . 13 -i 75.12 A (Pr=l) 

225648.42* 

305.6 

1.6 xlO - 4 

74 

0.163 

N 


269 , 18 - 268,19 A (vt—0) 

225648.42* 

261.7 

1.6 xlO ” 5 

106 



CH 3 OCHO 

21o,21-20i, 2 o A (ll7’ = l) 

226381.36* 

309.6 

2.8 xlO ” 5 

86 

0.277 

Y 


21i, 2 i-20i, 2 o A (v T =\) 

226382.72* 

309.6 

1.7 xlO ” 4 

86 




21 o. 2 i- 20 o. 2 o A (ur=l) 

226383.86* 

309.6 

1.7 xlO ” 4 

86 




21 i, 2 i- 20 o, 2 o A ((tr=l) 

226385.15* 

309.6 

2.8 xlO " 5 

86 



CH 3 OCHO 

21 o. 2 i- 20 i. 2 o E (yj=l) 

226433.26* 

308.9 

2.7 xlO ” 5 

86 

0.296 

N 


21 1 . 21 — 20 i, 2 o E (f 7 ’=l) 

226434.47* 

308.9 

1.7 xlO ” 4 

86 




21 o, 2 i- 20 o ,20 E (vj-1) 

226435.52* 

308.9 

1.7 xlO " 4 

86 




259,16-258,17 A (vt—0) 

226435.52* 

246.2 

1.6 xlO " 5 

102 




21 i, 21 - 20 o ,20 E (d7’ = l) 

226436.66* 

308.9 

2.7 xlO " 5 

86 



CH 3 OCHO 

2 O 2 . 19 —19 2 ,i8 A (vt—0) 

226718.69 

120.2 

1.7 xlO " 4 

82 

0.505 

Y 

CH 3 OCHO 

19 2 ,i7—18 2 ,i6 E (vt— 0) 

227019.55* 

116.6 

1.7 xlO ” 4 

78 

0.561 

N 


259 , 17 - 258,18 A (v T = 0) 

227021.13* 

246.2 

1.6 xlO " 5 

102 



CH 3 OCHO 

19 2 ,i 7 —1 82.16 A (vt-0) 

227028.12 

116.6 

1.7 xlO ” 4 

78 

0.436 

Y 

CH 3 OCHO 

2 O 4 . 17 -i 94 .i 6 A (t>r=l) 

242610.07 

321.7 

2.1 x 10" 4 

82 

0.092 

Y 

CH 3 OCHO 

37 7 ,3i-37 6 ,32 A (vt-0) 

242870.39* 

452.0 

1.8 xlO " 5 

150 

0.714 

N 


I 95 . 14 -I 85.13 E (vt-0) 

242871.57* 

130.5 

2.0 xlO ” 4 

78 



CH 3 OCHO 

I 95 . 14 -I 85.13 A (v T = 0) 

242896.02 

130.4 

2.0 xlO " 4 

78 

0.666 

Y 

CH 3 OCHO 

21 12 , 9 - 2111,10 E (v/ =0) 

316742.00* 

231.8 

3.3 xlO " 5 

86 

0.072 

Y 


21 i2,io-21 11,11 E (vt-0) 

316742.71* 

231.8 

3.3 xlO ” 5 

86 



CH 3 OCHO 

21 12 , 9 — 2111,10 A (vt-0) 

316776.74* 

231.8 

3.3 xlO " 5 

86 

0.069 

Y 


2112,10—21 11,11 A (vt-0) 

316776.74* 

231.8 

3.3 xlO ” 5 

86 



CH 3 OCHO 

26i3,i3-25i3,i2 E (t7y = 1) 

317177.16 

506.5 

3.6 xlO ” 4 

106 

0.043 

Y 

CH 3 OCHO 

98 , 2 - 87.1 A (ur=l) 

318009.06* 

256.8 

6.8 xlO ” 5 

38 

0.200 

N 


98 , 1 - 87,2 A (t>r=l) 

318009.06* 

256.8 

6.9 xlO " 5 

38 




13i2,i-13n,2 E (vt-0) 

318009.55* 

149.2 

1.3 xlO " 5 

54 




26i3,i3-25i3,i2 A Or=l) 

318012.17* 

506.2 

3.7 xlO " 4 

106 




26i3.i4-25i3.i3 A (itr=l) 

318012.17* 

506.2 

3.7 xlO " 4 

106 



CH 3 OCHO 

13 i 2 , 2 —13 11,3 E (vt-0) 

318016.90* 

149.2 

1.3 xlO " 5 

54 

0.155 

N 


98 , 1 - 87,1 E (vt-0) 

318017.37* 

69.0 

6.7 xlO " 5 

38 



CH 3 OCHO 

12i2,o—12 i 1,1 E (vt-0) 

318064.54* 

141.6 

7.0 xlO ” 6 

50 

0.204 

N 


98 , 2 - 87,1 A (vt-0) 

318065.26* 

69.0 

6.8 xlO ” 5 

38 




98 , 1 - 87,2 A (vr-0) 

318065.26* 

69.0 

6.8 xlO ” 5 

38 



CH 3 OCHO 

274 , 23 - 265,22 E ( Vt=1) 

318139.11* 

426.4 

3.4 xlO " 5 

110 

0.192 

N 


256 ,i 9 - 24 6 ,i 8 A (v r =l) 

318140.72* 

405.1 

4.6 x 10 - 4 

102 




2611 , 15 — 25 n,i 4 E (or=l) 

318145.25* 

474.5 

4.0 xlO " 4 

106 



CH 3 OCHO 

274 , 24 - 264,23 E (vt=1) 

318979.14 

417.8 

4.8 x 10 " 4 

110 

0.208 

Y 


Notes: The symbol * present after some frequency values indicates that the associated transition is blended with one or more 
transitions from the same species. a Y indicates that the line was considered in the rotational diagram analysis, while N indicates 

Article number, page 10 of 1 11 1 




A. Coutens et al.: Detection of glycolaldehyde towards the solar-type protostar NGC 1333 IRAS2A 


that it could not be used (for blending reasons). ^Blended with an unidentified species. ^Blended with CH 3 OCHO. ^Blended with 
CH 2 OHCHO. f Potentially blended with the gGg' conformer of ethylene glycol. 
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